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ABSTRACT 



Context. The pulse profiles of the transient Be/X-ray binary EXO 2030+375 show strong dependence on energy, as well as on its 
luminosity state, and are asymmetric in shape. 

Aims. We want to identify the emission components of the two magnetic poles in the pulsed emission to understand the geometry of 
the neutron star and its beam pattern. 

Methods. We utilize a pulse-profile decomposition method that enables us to find two symmetric pulse profiles from the magnetic 
poles of the neutron star. The symmetry characteristics of these single-pole pulse profiles give information about the position of the 
magnetic poles of the neutron star relative to its rotation axis. 

Results. We find a possible geometry for the neutron star in EXO 2030+375 through the decomposition of the pulse profiles, which 
suggests that one pole gets closer to the line of sight than the other and that, during the revolution of the neutron star, both poles 
disappear behind the horizon for a short period of time. A considerable fraction of the emission arises from a halo while the pole 
is facing the observer and from the accretion stream of the other pole while it is behind the neutron star, but the gravitational line 
bending makes the emission visible to us. 

Key words. Stars: neutron - X-rays: binaries 



1. Introduction 

EXO 2030+375 is an accreting X-ray pulsar with a pulsa- 
tion period of ~42 s, which wa s discovered with E XOSAT 
in 1985 during a giant outburst dParmar et alJ Il989bl) . A B0 
Ve star was found as its counterpart in follow-up observa- 
tions in the optical and infrared bands dJanot-Pacheco etail 
1988; iMotch & Janot-Pachecol fl987t ICoe et al.lll988l) . During 
the giant outburs t, EXO 2030+375 showed a spin-up of 
-P/P « 30 yr dParmar et all [l989b) a nd quasi-periodic os- 
cillations with a frequency of ~0.2 Hz dAngelini et al ] 119891) 
interpreted as caused by the formation of an accretion disk. 
Detailed analyses have shown that its rate of pulse-period 
change P, energy spe ctrum, and pulse pr ofil e are strongly 



lumin osity dependent ( Parmar et al . 1989 alhl: iRevno 



ds et alJ 



119931) . The orbital period is 46 days dWilson et al.1 120021 
and a normal outburst has been detected for nearly ev- 
ery periastron passage since 1991 dWilson et al.l 120051) . In 
2006, EXO 2030+375 u nderwent the first giant outburst since 
its di s covery in 1985 ( Corbet & Levind 120061 : iKrimm et alJ 
2006; iMcCollough et alJ 120061) . during which it reached a 
luminosity of Li-20kev ~ 1.2 x 10 38 erg s _1 
and again showed a strong spin-up. Rossi 



maximum luminosity 
dKlochkovetal.ll2008t) 
X-ray Timing Explorer (RXTE) monitored EXO 2030+375 ex - 
tensively during the 2006 giant outburst (Wil son et alJ 12008). 
The source was also observed by the IN TErnational Gamma 
Ray Astrop hysics Laboratory (JN TEGRAL. IWinkler et al.l2003l) 
and Swift dGehrels et al.l 120041). The spectra i n dicate a cy- 
clotron absorption line dKlochkov et al.l 120071 IWilson et al.l 



2008). Klochkov et al.1 d2008l) have shown that the spectrum of 
EXO 2030+375 changes with pulse phase, suggesting a fan 
beam geometry during the maximum, while towards the end of 
the giant outburst, it changes to a combination of a fan beam and 
a pencil beam. 

In X-ray pulsars, a neutron star accretes matter from a com- 
panion star via stellar wind or Roche lobe overflow. The ac- 
creted matter is channeled along the field lines of the strong 
magnetic field of the neutron star onto the magnetic poles. X- 
ray emission from the neutron star is produced in regions around 
the two magnetic poles. As the magnetic dipole axis is most 
likely inclined against the rotation axis of the neutron star, a 
distant observer sees pulsed emission. X-ray pulsars exhibit a 
wide variety of pulse shapes that differ from source to source. 
Generally, high -energy pulses have simpler shapes than low- 
energy pulses dWhiteetalJ 119831; iFrontera & Dalfiumd [19891: 
Bildst en et alJ 1 119971 and references therein). If one assumes an 
axially symmetric geometry for the two emission regions of the 
neutron star in a dipole configuration, the observed pulse profile 
should be symmetric. However, the observed pulse profiles typ- 
ically show an asymmetry. To explain the asymmetric shape of 
the total pulse profile, a distorted magnetic dipole field in which 
the two magneti c poles are not located opposite ea c h other have 
been discussed dParmar et al. 1989ah ILeahv 11991b [Riffert et al. 
119931: iBulik et al.lll995l) . iKraus et all d 19951) shows that, start- 
ing from the observed, asymmetric pulse profile, it is possible to 
disentangle the contribution of the two emission regions of the 
neutron star. Once the pulsed emission from each of the poles 
has been obtained, one can derive the geometry of the neutron 
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star. This again allows us to construct the beam pattern, i.e., 
the flux distribution from one emission reg i on. U sing this pulse- 
profile decomposition method. lKraus et all (11 996) have analyzed 
the pulse profiles of Cen X-3 and find indications of both pencil 
and fan beam. In the case of Her X-l, the resul ts of the pulse- 
profile decomposition by iBlum & Krausl (|2000) have not only 
shed light on the beam pattern of the magnetic poles, but have 
also confirmed that a warped and tilted accretion disk attenuates 
the emission from one pole of the neutron star. For A 0535+26, 
the reconstructed beam pattern suggests that the emission comes 
from a hollow colu mn plus a halo of scat tered radiation on the 
neutron star surface dCaballero et al.ll20Toh . 

In this paper we present the analysis of the energy-resolved 
pulse profiles of EXO 2030+375 utiliz ing the decomposition 
method developed by iRraus et aH dl995h . Section [2] gives an 
overview of the data used for our analysis and Sect.[3]describes 
the analysis and the results obtained with the pulse-profile 
decomposition method. The results are discussed in Sect. [4] 
Section [5] summarizes the possible geometry of the neutron star 
and the origin of the observed emission. 



2. Data 

EXO 2030+375 experienced a giant outburst in 2006, dur- 
ing which the source was monitored continuously by RXTE 
and was also observed by INTEGRAL. We have used 
the pulse profiles obtained with th e Joint European X-Ray 
Monitor (JE M-X. iLund et all [2003b and the imaging system 
IB IS/ISGRI (lUbertini et al.ll2003l) as presented in Figs. 2 and 8 
of lKlochkov et al.l(l2008h . 



2.1. RXTE observations 

For better statistics, we also used publicly available archival data 
from two observations with RXTE during the rise and the de- 
cay of the giant outburst when EXO 2030+375 showed about 
half of the maximum luminosity. The observations took place on 
June 28 and September 17, 2006 (observation IDs 91089-01-07- 
00 and 91089-01-19-01) when the luminosity of the source was 
£i-20kev = 7 x 10 37 erg s _1 and 6 x 10 37 erg s _1 , respectively. 
We started from the event files obtained with the Proportional 
Counter Array (PCA. Hahoda et al.ll 19961) . We used the event en- 
coded mode files in oder to have optimum binning. After fil- 
tering good time intervals and applying bitmasks, we created 
lightcurves with a time binning of 0.125 s in the following spec- 
tral bands: 2-9 keV, 9-14 keV, 14 - 20 keV, 20 - 30 keV, 
and 30 - 115 keV. After background subtraction, the lightcurves 
were corrected to solar barycenter and for orbitary motion of 
the binary. After folding the lightcurves with periods measured 
for each observation, we obtained pulse profiles with 128 phase 
bins. Phas e 0.0 was fixed t o agre e with the INTEGRAL pulse 
profiles of Klochk ov et al.l d2008). The folded lightcurves from 
the RXTE observations are shown in Fig.[T] 

During the two RXTE observations shown here, the lumi- 
nosity of EXO 2030+375 was comparable. As can be seen in 
Fig.Q] the pulse profiles before the maximum and after the max- 
imum of the giant outburst are very similar, corroborating that 
the s hape of the pulse pro files only depends on the luminosity 
state dParmaretal.lll989ah . 



3. The analysis 

3.1. The method 

A detailed description of the pulse-profile decomposition 
method, which is based on a back ward tracing of the emission, 
can be found in lKraus et al.l ([1995). All major steps of our anal- 
ysis and the criteria applied to obtain the best solution are de- 
scribed in the Appendix. 

The basic assumption of the method is that the magnetic 
dipole field of the neutron star is distorted in such a way that 
the two magnetic poles do not lie on a straight line through the 
center of the neutron star. Therefore, even though the emission 
from each pole is axisymmetric, the sum of the emission from 
both poles results in an asymmetric pulse profile. Using Fourier 
analysis, we model the observed asymmetric pulse profiles with 
two symmetric functions /12 to search for symmetry points O12 
in the pulse profiles and their offset A (see Appendix for de- 
tails). For each observation and energy range, the functions f\ 
and f2 correspond to the two single-pole pulse profiles that in 
total add up to the observed asymmetric pulse profile. Each sym- 
metry point corresponds to the pulse phase where the respective 
pole is either closest to or most distant from the observer's line 
of sight. From the two symmetry points and functions, we then 
derive the location of the emission regions and the beam pattern. 

3.2. Decompositions 

We have a total of 26 pulse profiles from INTEGRAL and RXTE 
observations. To perform pulse-profile decomposition, the max- 
imum in each of the 26 pulse profiles is normalized to unity. 
Thereafter each pulse profile is modeled with two symmetric 
functions f\ and based on Fourier analysis. After applying 
criteria 1 (positive flux, see Appendix lA.il ) and 2 (no ripples), 
we obtain a large number of possible values for the parameters 
<t>i and A. However, after combining the results for all pulse pro- 
files (criterion 3), only two interesting solution regions remain in 
the parameter space of <t>i-A. We call these solutions 1 and 2 and 
perform further analysis with these two possible solutions. 

For each observation and energy band, the functions f\ and 
f% correspond to the two single-pole pulse profiles that in total 
add up to the observed asymmetric pulse profile. The Fourier 
analysis finds more than one possible set of 3>i and A within a 
small region for one solution. For each total pulse profile, we 
have to look at the different sets of the single-pole pulse profiles 
to decide which one is consistent with the single-pole pulse pro- 
files at other energies. The pulse profiles of the different obser- 
vations are studied separately, because one should see a correla- 
tion between the different energy bands of one luminosity state, 
but not necessarily between two different observations. Figure 
[2] shows some of the selected single-pole pulse profiles and the 
derived beam patterns for each pole for solution 1 . Within one 
observation, one can see an energy-dependent evolution of the 
single-pole pulse profiles. The values of <l>i and A for solutions 
1 and 2 are (65° < <J>i < 75°, 63° < A < 70°) and (70° < 4>i < 
80°, 81° < A < 88°), respectively. The parameters for the two 
solutions are listed in Table Q] 

3.3. Overlaying beam patterns 

As the neutron star rotates, the angle between the axis through 
one magnetic pole and the line of sight 9 changes with phase, i.e., 
with rotation angle O. The decomposition has provided us with 
beam patterns as seen by the distant observer for each emission 
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Fig. 1. Pulse profiles of EXO 2030+375 during the rise (left) and the decay (right) of the giant outburst in 2006 observed with 
RAT/:. 

Solution 1 




Phase (*/2it) 




Fig. 2. A selection of single-pole pulse profiles (left) and corresponding beam patterns (right) obtained by the decomposition 
method for the O^A set for solution 1 for the energy bands 3-9 keV, 9-14 keV, 14 - 20 keV, 23 - 26 keV, and 39 - 45 keV of 
the INTEGRAL observation near maximum. The solid line is used for emission from the first pole, dashed line for the second pole. 
The vertical dash-dotted lines indicate the ranges for the symmetry point fl>i, while the dotted lines show the ranges for <t>2 obtained 
from pulse profiles at different energies and observations. 



region as functions of the phase. Now we compare the two beam 
patterns derived from the two single-pole pulse profiles of each 
observed pulse profile and search for a range in cos(<l> - <$>i), and 
thus 6, in which the two beam patterns seem to show the same 
emission (see Appendix IA.2t . We try to overlay the two beam 
patterns by using the relation: 

cos(d> - $0 = a + b cos(0 - 2 ), b > 0. (1) 



For solution 1 we are not able to find an overlap of the single- 
pole beam patterns. It is more likely that, in this case, the geome- 
try only allows the observer to see two different parts of the total 
beam pattern. To assemble the two parts of the beam pattern, we 
use a = -2.1 and b — 1.0 (see Fig. IA.21 upper panel). For so- 
lution 2, we find an overlap of the single-pole beam patterns for 
a = -0.1, b = 1.0, although the beam patterns do not seem to 
match perfectly. 
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Fig. 3. Asymptotic beam patterns of the plus and minus solutions of solution 1 seen by the distant observer plotted over the angle 
9 between the first magnetic pole and the line of sight. The left panel shows the beam patterns of solution 1 of the observation near 
maximum as shown in Figure [2] The right panel shows a selection of additional data from the observation after maximum with 
similar energy bands. 



Table 1. Parameters obtained from the decomposition for solu- 
tions 1 and 2 





A 


a 


b 


0! 1 


2 ' 


S l 


1 65° -75° 


63° - 70° 


-2.1 


1.0 


39° 


141° 


40° 


2 70° - 80° 


81° -88° 


-0.1 


1.0 


87° 


93° 


85° 



1 Assuming &o = 50°. 



A total beam pattern can be reconstructed from the beam pat- 
te rns calculated from the single-pole pulse profiles. As shown 
in Kraus et al .1 d!995h . there is an ambiguity in the relation be- 
tween O and 9, as each single-pole pulse profile has two sym- 
metry points at O,- and O,- + n. Therefore, the solutions cannot 
tell us which ends of the sections of the beam patterns belong 
to, e.g., the lower values of 9. For each set of Oi and A we ob- 
tain two possible solutions (called plus and minus) for the total 
beam pattern. The decomposition method cannot tell us which 
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Fig. 4. Relation between the observing angle ©o and the angles 
©i (dotted line), ©2 (dashed line), and 5 (range of possible val- 
ues) for solution 1 . 



one is the real solution. We have to take results from other mea- 
surements of the source into consideration, e.g., luminosities or 
spectra that will give insight into the emission processes, to de- 
cide which of the two is the real solution. The reconstructed 
asymptotic beam patterns for solution 1 as seen by the distant 
observer are shown in Fig. [3] Here the beam patterns are plotted 
over the angle 9 between the normal at the first magnetic pole 
and the line of sight. 

That the beam patterns of the single-pole pulse profiles do 
not match in the overlap region of solution 2 might indicate that 
the emission from the two magnetic poles are not identical and 
cannot be described with one local beam pattern. Therefore, we 
take the mean of the two beam patterns and model a total aver- 
aged beam pattern, which we use to reconstruct the visible to- 
tal pulse profile, i.e., assuming equal local emission pattern for 
the two poles. The reconstructed pulse profiles (Fig. IA.3b show 
significant deviations from the observed profiles. In addition, 
this solution yields a very extreme geometry with the two mag- 
netic poles located near the equator of the neutron star, forming 
an angle of ~90° between each other (see Table [TJ. Although, 
in principle, we cannot rule such a strongly distorted geome- 
try out, especially not for young neutron stars like those ex- 
pected in Be/X-ray binary systems, as they might have experi- 
enced some anisotropic conditions while their birth, this rather 
unlikely geometry also suggests that solution 2 is not appropri- 
ate for EXO 2030+375. Therefore, in the following, we focus on 
the discussion of solution 1. 



3.4. Geometry of the neutron star 

To derive the exact geometry of the neutron star, i.e., to de- 
termine the polar angles of the magnetic poles ©i and ©2 as 
well as the offset angle 6, we need to know the inclination an- 
gle of the rotation axis of the neutron star. However, this angle 
is not known for most of the neutron stars, in particular not for 
EXO 2030+375. Figure [4] shows the dependence of the angles 
©i, ©2, and S on the angle ©o between the rotation axis of the 
neutron star and the line of sight of the observer. To convert the 
phase parameter cos(0 - <&i) into the angle 9 between the first 
magnetic pole and the observer, we have to assume an inclina- 
tion angle ©o with respect to the rotation axis. Here, we use ©o = 




Fig. 5. Decomposed single-pole pulse profiles of the data taken 
near the maximum of the giant outburst for the lowest energy 
band (same as the upper left diagram of the upper [solution 1] 
panels in Fig.|2| with proposed emission components derived 
from the beam patterns of the minus solution (Fig. [3] lower 
panel). We suggest that the emission from the first pole con- 
sists of emission from the halo (diagonally striped component) 
and from the accretion column (vertically striped). The domi- 
nant emission from the second pole seems to have its origin in 
the scattered photons in the upper accretion stream, which is vis- 
ible while the second pole is behind the horizon of the neutron 
star due to gravitational light bending, while there is also a con- 
tribution from the accretion column (vertically striped). In this 
solution, the first pole gets closer to the line of sight than the 
second. 

50°, corresponding to ©i = 39°, © 2 = 141°(see Fig.|4]and Table 

4. Discussion 

In this section we want to further examine the results of the de- 
composition method and present the possible geometry of the 
neutron star. 



4.1. Disentangling the emission components 

For solution 1, there is no overlap between the two beam pat- 
terns obtained from the single-pole pulse profiles as shown in 
Sect. 13. 31 It means that only a part of the emission is seen from 
each pole during the revolution of the neutron star. By putting the 
two visible parts together, we obtain the total beam pattern of the 
emission around one magnetic pole. In doing so, we make the 
assumption that the two magnetic poles have the same emission 
pattern. The plus solution (Fig. [3] upper panel) can be described 
as a composition of a forward directed emission (towards 0°) that 
is more pronounced at higher energies and an extended, rela- 
tively soft emission (~80° - 180°). The harder, forward-directed 
emission indicates a pencil beam. The relative flux of the softer 
emission component is higher in the data from the observation 
performed near the maximum compared to the one during the de- 
cay and can be interpreted as a fan beam, in agreement with the 
largely accepted picture that an optical ly thick accretio n column 
is formed during the giant outburst (IWhite et alJI 1983L and ref- 
erences therein). The minus solution has a soft emission at 9 $ 
60°and a harder emission at 9 ^ 120° that increases for larger 
9. As newest calculations by Kraus et al. (in prep.) have shown, 
reprocessing of photons in the upper accretion stream creates a 
significant emission component that dominates the beam pattern 
at higher energies and can be observed while the pole, hence the 
accretion column, is on the other side of the neutron star ('anti- 
pencil'). In addition, the emission from a halo that is formed by 
scattered photons at the bottom of the accretion column domi- 
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nates the beam pattern at lower energies and at lower 8 and is 
stronger in the data near the maximum of the giant outburst. 

Figure [5] shows how the emission from the two poles con- 
tributes to the pulse profiles for the decompositions of the low 
energy band data taken at the maximum of the giant outburst 
in the case of the minus solution. Let us assume that the first 
magnetic pole with the polar angle @i gets closer to the line of 
sight than the second magnetic pole, i.e., ©o - ©i < ©2 - ©o- At 
phase 0.0, both poles are right behind the horizon of the neutron 
star: the first pole is going to reappear, the second pole is turn- 
ing farther away from the observer. At phase -0.1 the first pole 
becomes visible, and from then on the emission from the halo 
of the first pole makes the largest contribution until the pole dis- 
appears behind the horizon at ~ 0.6. When the first pole, which 
is closer to the observer's line of sight, is right at the horizon of 
the neutron star and its accretion column is seen from the side, 
a minimum is likely to be observed in the pulse profile (see also 
Sect. l4.4l >. At the major maximum at phase 0.15 - 0.25, the sec- 
ond magnetic pole is behind the horizon of the neutron star and 
the scattered and gravitationally bent photons from the upper ac- 
cretion stream cause the pronounced increase in flux. The second 
pole that is rotating on a circle farther away from the line of sight 
than the first pole comes back to the front side at about phase 
0.45 and is closest to the observer at about phase 0.7. However, 
since the line of sight is closer in latitude to the first pole than to 
the second pole, the observer never gets as close to the surface 
normal of the second pole as to that of the first pole. The main 
emission seen from the second pole while it is on this side of the 
horizon comes directly from the accretion column. 

iParmar etal . (1989a) have modeled the luminosity depen- 
dent pulse profiles from the first observed giant outburst in 1985 
by assuming a fan beam and a pencil b eam component for th e 
two magnetic poles based on a model bv lWang & Welted (Il98ll) . 
They obtain a fit for all ten pulse profiles for different luminosi- 
ties ranging from 0.1 - 10.0 xlO 37 erg s with some residuals 
and find that the fan beam mainly produces the major peak at 
phase 0.2, whereas the peak at phase 0.95 can be ascribed to a 
pencil beam. The way they have chosen phase 0.0 is different 
than in our w ork, resulting in a sh ift by about 0.4. The best-fit 
parameters of Par mar et al.l (Il989al) correspond to ©i = 70°, ©2 
= 1 10°, and 6 = 70°, assuming ©o = 25°. This result is consistent 
with our solution 1 for the case of © = 25° (see Fig.©. 

4.2. Comparison to phase-resolved spectral analysis 

iKlochkov etafl d2008l) analyzed the same INTEGRAL data as 
used in this work and performed pulse phase resolved spec- 
troscopy. For the data of the maximum, they find that the spec- 
trum is harder at the main peak (phase 0.2). They conclude that 
the observer might be seeing Compton scattered photons from 
the optically thick accretion column. In the data taken at the end 
of the decay, a new peak appears at about phase 0.95. At this 
point, the spectrum is again harder, leading to the conclusion that 
the line of sight of the observer is closest to the magnetic field 
lines. This picture is in good agreement with the solution 1 of 
the decomposition, in which the emission of the upper accretion 
stream of the second pole, which is visible owing to gravitational 
light bending, has its maximum at phase ~0.2 while the halo of 
the first pole also has substantial emission, thus contributing sig- 
nificantly to the main peak. At phase ~0.95 the scattered emis- 
sion from the upper accretion stream of the second pole starts 
to increase, while the emission from the accretion column of the 
first pole also contributes to the observed emission. 



4.3. Intrinsic beam patterns 

The observed asymptotic beam pattern differs from the intrinsic 
local beam pattern of the emission region because of relativistic 
light deflection near the neutron star. Making assumptions on 
the radius r and mass m of the neutron star, the asymptotic angle 
6 can be transformed into the intrinsic angle i? relative to the 
surface normal at the first magnetic pole to describe the local 
emission pattern. We use the canonical values r — 10 km and 
m = 1.4 Mq, thus the ratio between r and the Schwarzschild 
radius rs is r/rs = 2.4. The intrinsic beam patterns of the minus 
solution are shown in Fig. [6] in polar diagrams. As can be seen 
in these polar diagrams, the range for the visible angle is ~10° - 
105° if ©o = 50° is assumed. 

4.4. Possible origin of the emission 

In X-ray pulsars, a strong magnetic field funnels matter accreted 
from the companion star onto the polar caps. The infalling par- 
ticles deposit their energy in the atmosphere of the neutron star. 
Heating due to accretion is balanced by radiative cooling through 
bremsstrahlung emission, Compton scattering, and the so-called 
cyclotron radiation after collisional excitation of electrons. In 
low-luminosity states, the accreted matter is decelerated in the 
atmosphere of the neutron star by Coulomb scattering. Emission 
can be seen as pencil beams from the hot spots at polar caps 
dWang & Welterlfl98ll) . In high-luminosity states, however, the 
accretion rate is higher, and an accretion column is believed to 
form. Plasma is decelerated by radiation pressure in the col- 
umn, and a radiative shock forms above th e neutron star surface 
(lDavidsonlll973t iBasko & Sunvaevlll976l) . Above the radiative 
shock, plasma is in free fall, while below the discontinuity, there 
is a region of nearly stagnant plasma from which photons escape 
from the sides of the column in a fan beam. 

In add ition, a luminous halo might f orm around the accretion 
colum n ([Davidson & Ostrikei] 119731; iLvubarskii & Svunvaevl 
1 19881) and radiation from the polar cap and the halo can also 
be scattered in the upper a ccretion stream dSoffel et"al] 1985: 
iBrainerd & Meszaroslll991l) . As the accretion stream is delim- 
ited by magnetic field lines of the neutron star, it opens up wide 
far from the neutron star. Therefore, the emission from the upper 
part of the accretion stream can dominate the observed flux and 
can also screen the polar caps and the halo. Due to relativistic 
light deflection, emi ssion from one pole can be deflected to the 
antipodal direction dPechenick et all Il983t iRiffert & Meszaroll 
1988: lLeahv & Lill 1 995t lKrausll200 lh . All these effects can mod- 
ify the local beam pattern a nd thus have an effect on the pulse 
profile of the X-ray pulsar. IWang & Welter! d!981l) have mod- 
eled the pulse profiles of X-ray pulsars by assuming hot spots 
and fan beams, although EXO 2030+375 is not included in their 
list of selected X-ray pulsars. They show that, in general, sharp 
minima can be seen in the pulse profiles at the moment when the 
fan beam rotates behind the horizon of the neutron star and when 
it reappears on the side of the neutron star facing the observer. 
This is exactly what seems to be happening in EXO 2030+375 
at phases 0. 1 and 0.6 as the decomposition method has shown. 

Kraus et al. (120031 in prep.) have modeled beam patterns and 
pulse profiles for medium-luminosity X-ray pulsars assuming an 
accretion column with energy-dependent local beaming of radi- 
ation, a luminous halo formed by illumination of the neutron star 
surface, and magnetic scattering in the upper accretion stream. 
They show that the emission from the upper accretion stream can 
dominate the local beam pattern for 9 > 120°. Halo emission has 
a maximum at 6 « 30° - 60°, with significant contribution to 
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Fig. 6. Polar diagrams of the intrinsic beam pattern for the minus solution of solution 1 for the same observations as shown in 
Fig. |3] for the softest and the hardest bands for an assumed ©o = 50°. Dotted lines are drawn to clearly show the range in which the 
emission is visible for the observer. 



the beam pattern for photon energies below ~5 keV. If one as- 
sumes isotropic emission from the accretion column below the 
discontinuity, the halo emission becomes comparable or negli- 
gible relative to the column emission at ^10 keV. However, if 
the radiation below the shock is beamed downwards, the halo 
remains more luminous than the accretion column and the upper 
stream up to ~30 keV. The decomposition has shown that the 
beam patterns show dominant soft emission for 6 < 60°, which 
can be interpreted as halo emission, while the emission 6 > 150° 
most likely arises from the upper accretion stream (see Fig. [3] 
lower panel). 

5. Summary 

We performed pulse-profile decomposition with data of 
EXO 2030+375 taken during the giant outburst of 2006 by 
RXTE and INTEGRAL. This is the fourth source after Cen X-3, 
HerX-1, and A 0535+26 to which this method has been applied. 
Each of the asymmetric pulse profiles of EXO 2030+375 at dif- 
ferent luminosities in various energy bands are decomposed in 
two symmetric pulse profiles that account for emission from the 
two emission regions of the system. 

We find that the magnetic field of the neutron star is mod- 
erately distorted. The observer sees a part of the emission from 
each of the two emission regions, but these parts do not over- 
lap. We suggest that the main peak at phase 0.2 in the observed 
pulse profiles can be attributed to harder emission from the up- 
per accretion stream of the second pole, which can be observed 
while the pole is on the other side of the neutron star because of 
relativistic light bending around the neutron star ('anti-pencil'). 
However, the main peak also has a considerable contribution 
from the halo emission of the first pole, which is closer to the 
line of sight than the second. Between phases ~0.95 and ~0.1, 
both poles are located on the rear side of the neutron star. The 
sharp minima seen at phases 0.0 and 0.6 are caused when the first 
pole is about to re- and disappear at the horizon of the neutron 
star, and its accretion column is seen from the side. 

Our analysis has disentangled the emission components of 
the neutron star, which in total lead to the observed asym- 
metric, energy-, and luminosity-dependent pulse profiles of 



EXO 2030+375. It will allow us to perform detailed analyses 
of, e.g., the pulse phase-resolved spectra with reliable interpreta- 
tion of the differences in the spectral parameters, hence shedding 
light on the physical processes in the system responsible for the 
observed emission. 
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Appendix A: The decomposition method 

The pulse-profile decomposition m ethod has been developed and 
first presented bv iRraus et al.l (Il995l) . Here we quickly summa- 
rize the method and the steps in the analysis. 

A. 1 . Decomposition into two single-pole pulse profiles 

Let 9 be the angle between the direction of the line of sight and 
the axis through one magnetic pole. Figure lATl shows the con- 
figuration of the neutron star, showing the positions of the mag- 
netic poles with respect to the rotation axis. The polar angles of 
the two magnetic poles are called @i and ©2. As the neutron star 
rotates, 9 changes with rotation angle <t>. For the emission of each 
pole, there are two symmetry points at i>, and <J>, + n (i = 1,2). 
For each pole, the relation between 9, ©,, and <J>, can be obtained 
using the spherical triangle in Fig. lA.ll (left): 

cos 9 = cos ©0 cos ©, + sin ©0 sin ©, cos(<J> - <£,). (A.l) 

We assume that there is an offset 5 between the second magnetic 
pole and the antipodal position of the first magnetic pole, there- 
fore, there will be a phase offset for the rotation angles <E>, from 
a symmetric configuration which can be defined as 

A = 7r - (<J>i - <D 2 ). (A.2) 

For an ideal dipole field ©i + ©2 = n, A = 0, and (5 = 0. With 
the polar angles 0j, ©2, and the offset A in the rotation angles, 
we have a complete set of parameters to describe the geometry 
of the neutron star. 

To find the contributions from each magnetic pole, we per- 
form Fourier analysis of the observed total pulse profiles. We 
model it as a sum of two symmetric functions f\ and fz and 
search for the values for their symmetry points <t>i and <D 2 , re- 
spectively. In principle we are able to find a set of f\ and fa for 
any chosen <l>i and <J>2- However, since we deal with functions 
that describe astronomically observed emission, the following 
criteria need to be fulfilled: 

1 . Positive flux: The symmetric functions f\ and / 2 must not 
have negative values because they model the flux of an as- 
tronomical object. 

2. No ripples: The functions f\ and should show no small- 
scale features that cancel out in the sum. As the two func- 
tions correspond to pulse profiles of single-poles that emit 



independently, they ought not to have features that match ex- 
actly. Also, the single-pole pulse profiles are supposedly not 
more complicated than the total pulse profile. 
3. Same geometry: We have pulse profiles from one source in 
different energy bands and, ideally, from more than one ob- 
servation. The observed pulse profiles are energy and lumi- 
nosity dependent in most cases. Since the emission arises 
from only one object, the symmetry points must be the same 
for the decompositions of all available data. 

The actual search for reasonable decompositions is performed 
using the parameters <J>i and A. The second symmetry point <1>2 
can then be obtained from Eq. IA.2l The decompositions of the 
pulse profiles of different energy bands and different observa- 
tions will not yield identical values for <J>i and A. Therefore we 
have to define interesting regions in the parameter space of <t>i 
and A and look at the single-pole pulse profile for similar val- 
ues of <t>i and A for each energy band and observation to decide 
if they can be declared as one consistent set of emission. This 
means that the single-pole pulse profiles have to have similar 
values of <J>i and A and ought to be similar from one energy band 
to the other in one observation. Once one has found such a set of 
decomposition for all energy bands and observations for similar 
<l>i and A, one can compute the corresponding beam patterns for 
each pole. 

A.2. Reconstruction of the beam pattern 

Two beam patterns are obtained from single-pole pulse profiles 
as functions of the rotation angles <D (for one single-pole pulse 
profile) and (for the other single-pole pulse profile). During 
one revolution of the neutron star, the angle 9 between the first 
magnetic pole and the line of sight changes with the phase. There 
may be an interval during one phase, i.e., a range of the angle 9, 
in which the observer sees emission from both poles. The same is 
true for the second magnetic pole. The range covered by the an- 
gle 9 is in general different for each pole. There may, however, be 
a certain subrange of values of 9 that occur for each of the poles. 
Then, at some phase <£, the observer looks onto the first pole at 
angle 9 and at some different phase O, the observer looks onto 
the second pole at the same angle 9. A simple example would be 
two antipodal poles that pass through the line of sight (9 = 0) at 
<J> = and at O = 0.5, respectively. If, in addition, the two emis- 
sion regions of the neutron star have the same beam pattern, the 
emission seen at 9 from the first pole (phase <D) is the same as the 
emission seen at 9 from the second pole (phase €>). This means 
that the visible beam patterns of the two poles must have iden- 
tical parts at different pulse phases. If we find such parts of the 
beam patterns of the two single-pole pulse profiles we can over- 
lay them and get a relation between cos(<J>-<Di) and cos(4>-<J>2)- 
From this relation, we can derive the positions of the magnetic 
poles ©i and ©2 as functions of the direction of the observation 
©0. The angle ©0 needs to be determined independently in other 
studies of the source. 

List of Objects 

'EXO 2030+375' on paged 




Fig. A.l. Intrinsic geometry of the neutron star. Figures taken from lKraus et al J d 1995b 
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Fig. A.2. Beam patterns from the two emission regions seen by the observer plotted over cos(cD-cDi). The left panel shows the beam 
patterns of solution 1 of the observation near maximum as shown in Figure [2] The right panel shows a selection of additional data 
from the observation after maximum with similar energy bands. For solution 1 (upper panel), In this case no overlap of the beam 
patterns for the two emission regions is found. For solution 2 (lower panel), the beam patterns seem to overlap over a wide range, 
but show small differences. For the parameters a, b, see Section [331 
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Fig. A.3. Left: Beam patterns from the two emission regions plotted over cos(<t> - <Di) for the energy band 26 - 30 keV of the 
INTEGRAL observation near maximum for solution 2. For this solution, the beam patterns seem to overlap over a wide range, but 
show small differences between each other. Right: Observed total pulse profile (data points with errors) with total pulse profiles 
reconstructed from the not perfectly matching beam patterns of solution 2 for the two emission regions (dotted) and averaged beam 
patterns assuming equal local emission pattern for the two regions (dashed). 



